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(SI ' The operator product expansion (OPE) in 4d (super) conformal field theory is of broad 

interest, for both formal and phenomeno logical applications. In this paper, we use con- 
formal perturbation theory to study the OPE of nearly-free fields coupled to SCFTs. 
Under fairly general assumptions, we show that the OPE of a chiral operator of dimen- 
sion A = 1 + e with its complex conjugate always contains an operator of dimension less 
than 2A. Our bounds apply to Banks-Zaks fixed points and their generalizations, as we 
^ ! illustrate using several examples. 
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1. Introduction 

Conformal field theories (CFTs) and superconformal field theories (SCFTs) play a 
valuable role in our understanding of quantum field theories. Many interacting theories 
are believed to flow to (S)CFTs, while many more can be understood as small deformations 
away from conformality. Furthermore, through the AdS/CFT correspondence, CFTs are 
a useful tool for understanding quantum gravity. 

Combined with unitarity, the enhanced symmetry of these theories has long been 
known to impose significant constraints on the spectrum of operator dimensions |I],£2J. 
These "unitarity bounds" have played a valuable role in our understanding of the phases 
of supersymmetric gauge theories @ . In addition, a-maximization [|] has made it possible 
to compute the dimensions of chiral operators in proposed IR SCFTs, and unitarity bounds 
can be used to test the consistency of such proposals. 

In recent years, there has been significant progress in using conformal invar iance to 
constrain the spectrum of primaries that can appear in the operator product expansion 
(OPE) of scalar operators. Specifically, given a CFT with a scalar primary operator X of 
dimension Ax, the OPE takes the form 

X\x)X{V) = -±- + £ ^^ + . . . , (1.1) 

i 

where the q are the OPE coefficients, Aj are the dimensions of the scalar primary operators 
Si and . . . signifies higher spin operators and conformal descendants. Applying the OPE 
to the four point function of X, [H-H] were able to place an upper bound on 



Smin = A m j n — 2Ax (1.2) 

where A m i n denotes the smallest Aj that appears in the OPE. Subsequent work has ex- 
tended these results, including to supersymmetric theories P~P1- 

In the context of supersymmetry, there has been much interest in theories with large 
anomalous dimensions, in part because of their many applications to model building, see 
e.g. [12-2SJ (for non-SUSY motivation, see e.g. |5,3D[ and references therein). For example, 



one can try to resolve the fi/Bfj, problem by using an SCFT with S m i n > [pT2| - p~5j] . We 
often imagine generating \x and B/i at the SUSY breaking scale through the F-term of a 
chiral operator X via the effective lagrangian 
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Electroweak symmetry breaking requires cg^ ~ c 2 at the scale SUSY breaking. Unfortu- 
nately, models almost invariably give rise to CBp. ^> ci at the scale where they are generated 
- this is known as the n/Bfj, problem (see [RT] for a review and original references). How- 



ever, if the theory flows near a SCFT in between these scales, the u X^X n is replaced by 
the operators in the OPE of AU and X. If alffl such operators have dimension Aj > 2Ax, 
then RG flow will suppress cg M down to acceptable levels in the IR. 

Although the bounds in P^TTf allow for positive 5 m i n , to date there is no existence 



proof showing that this is possible. In []33[], SCFTs with weakly coupled AdS duals were 
studied, and it was shown that Srnin could have either sign as far as the 5d effective 
gravity theory was concerned. If these models could be UV completed into full string 
compactifications, then they would furnish existence proofs of positive 5 m i n . Whether this 
is possible remains an interesting open question. 

In this paper, we will take a complementary approach, and investigate the quantity 
S m i n in the context of SCFTs in which X is a chiral primary operator with dimension 
Ax = 1 + e, with e<C 1. When Ax = 1, the unitarity bounds require that X is a free 
field and its OPE is trivially determined. One might suspect that for Ax = 1 + e, the 
bounds on operators in the OPE are accessible perturbatively in e. We will develop a 
general formalism for computing 5 m i n based on conformal perturbation theory, with the 
assumption that X acquires its anomalous dimension solely through a coupling in the 
superpotential to a chiral operator O of dimension A© = 2 — e. Under such circumstances, 
the theory will flow to a fixed point where X has the required dimension, and we will prove 
that the OPE of X and X^ contains an operator of dimension at most A = 2Ax + |C(e 2 )|. 
We will further show that S m [ n ~ —e 2 /i>i < when the original SCFT contains real 
operators Li of dimension Aj = 2 + z/j with ^ < 1. As an application of our formalism, 
we will prove that in fully perturbative examples, such as generalized Banks-Zaks fixed 
points, one always finds d m in < 0. 



One may weaken this requirement to allow conserved current multiplets J (i.e. real multiplets 
of dimension two) to appear in the OPE |3^| . This is viable because D 2 J = 0, by definition, and 
therefore does not contribute to Bfi. However, one cannot spontaneously break SUSY in an 
SCFT. Therefore, the SCFT must be deformed by relevant operators in order to break conformal 
invariance. For this to be a valid solution to the \i — Bfi problem, one must show that the 
deformation can be chosen such that the current remains conserved at all energies. We will ignore 
this possibility, as it requires knowledge of the field theory that cannot be determined from the 
behavior near the fixed point alone. 



While these assumptions may seem restrictive, it is plausible that this setup actually 
applies to all examples of SCFTs with a chiral operator with dimension close to one. In 
the limit e — > 0, the operator X must become free and will decouple from any other fields. 
Furthermore, X must be gauge invariant. Provided there is some Lagrangian description 
of X, this suggests the above description should hold in some duality frame. 

The organization of the paper is as follows: In Section 2, we will explain how the 
bounds on scalar operators in the OPE of chiral primaries in an SCFT can be under- 
stood from conformal perturbation theory. The presentation in this section will be mostly 
general, making only a few simplifying assumptions about the form of the SCFT. In Sec- 
tion 3, we will compare these results with concrete examples derived from Banks-Zaks-like 
SCFTs. The results in this section agree with the results of conformal perturbation theory 
but are derived in a more conventional language. We will summarize our results and men- 
tion some interesting open problems in Section 4. Appendix A contains the details of the 
conformal perturbation theory calculations. Appendix B addresses complications arising 
from the presence of global symmetries. Appendix C contains further details regarding the 
matching between the results of Section 2 and 3. 

2. Scaling Behavior of Nearly Free Fields 

2.1. Setup 

Consider a 4d interacting SCFT V\ which contains a chiral operator O of dimension 
A = 2 - e. The OPE of O and 0+ is given by: 

C(x) f C(0) = \x\~ 2A ° + J2 \x\ Ai ~ 2Ao c,L, + ... , (2.1) 

i 

where Lj are real scalar multiplets with dimension A* = 2 + z/j, Vi > 0, and . . . denotes 
descendants and operators with higher spin. All primary scaling operators in this paper 
will be "CFT-canonically normalized" as in ( |2.1| ), i.e. with the coefficient of the leading 
singularity set to one. 

Now consider adding to V\ a free chiral superfield X, and the deformation 

c = c Vl + -^ fd 4 ex^x+ ( fd 2 e— A e ox + h. c . j , (2.2) 

where we have introduced the Wilsonian cutoff scale A in order to make A dimensionless. 
Note that X is not a canonically normalized field; this is to obey the more natural CFT 



normalization described in the previous paragraph. When < e < 1, this theory will 
flow to a new fixed point, V2, with A ~ C^i/e). At the new fixed point, we would like to 
know what is the lowest dimension scalar operator that appears in the OPE of X and X^ . 
Candidate operators include X^X and the L^. In general, these are not scaling operators 
at the new fixed point V2, but will instead mix with one another. As we will see in the 
following subsections, when the dimensions of Li are approximately 2, this mixing is the 
dominant contribution to £ m in- 

We will focus on the case where conserved currents are absent in the X^(x)X(0) OPE. 
In the free theory, X^X is the conserved current associated with the symmetry U(l)x- 
This current is broken by ( |2.2| ) but, in general, new conserved currents could arise from a 
combination of U(l) currents broken by the superpotential interaction, i.e. a combination 
of U(l)x and C/(l)'s from V\ under which O is charged. The unitarity bounds ensure 
that any such current would have dimension two, and therefore would trivially satisfy 
A min < 2 A x for Ax > 1. 

To simplify the discussion in this section, we will assume that O is a singlet under all 
the global symmetries of "Pi, i.e. the Vi are strictly positive in ( |2.1| ). This ensures that the 
deformation ( |2.2| ) breaks the U(l)x symmetry completely, and then absence of conserved 
currents is trivially guaranteed. In appendix B, we will generalize to the (considerably 
more complicated) case of charged O. With some mild assumptions, all the results of this 
section carry over to charged O, provided that conserved currents remain absent in the 
X*X OPE. 

2.2. RG flow 

Under the assumption that e C 1, we can understand the RG flow induced by (|2.2| ) 
using conformal perturbation theory. Without loss of generality, A will be real, as its phase 
can be removed by a field redefinition of X . As above, A will denote the sliding Wilsonian 
cutoff scale. 

Once we introduce the deformation by a single operator, as in ( |2.2j ), other operators 
will also contribute to the RG flow at higher orders in A. Therefore, the RG is actually 
described by the flow in a larger space of couplings. We will include these additional 
operators in the Lagrangian as 

C = C Vl + fd 4 6 (^Z x X^X + ^y i K-^L i + ..\ + ff d 2 6 ^A e OX + h.c^\ . 

(2.3) 



Here Zx is the wavefunction renormalization of X and yi are coupling constants. The . . . 
denote operators that that do not appear in the 0*0 OPE. The counterterms for these 
can only arise at higher orders in A, and so they will not contribute to our analysis in 
this paper. We have also dropped operators of the form X^XLi as they contribute at the 
same order as Li of dimension A« > 4 and will be negligible in the regimes of interest. 
We have included explicit A dependence to make all the couplings dimensionless. By a 
non- renormalization theorem ||34j| , we have chosen a holomorphic scheme such that no 
additional operators appear in the superpotential. 

If we work in a holomorphic scheme, there is no renormalization of A and its beta func- 
tion (5\ = d Yo R \ is given purely by the classical contribution, (3\ = — eA. In addition, there 

( \ \ 2lx ( x > yi } 
is wavefunction normalization, namely Zx — ( a ) = l+27x(A, yi) log(Ao/A) + . . . 

(where Aq is some arbitrary reference scale). In this scheme, as we flow to the IR (A — > 0), 

A and Zx will run to infinity. To see that we reach a fixed point, we redefine X and A to 

make the former (nearly) canonical 

x _> x phys = z]/ 2 X, A -> X phys = Z~ 1/2 X . (2.4) 

After the redefinition, the beta function of A (henceforth dropping the "phys" superscripts) 
becomes 

& = -eA + A 7 x(A,7/i) . (2.5) 

Regardless of the form of 7x(A, yi), we see that the fixed point will arise at 7x(A*, y^*) = e. 
As expected, the dimension of X at V2 is Ax = 1 + 7x(A*, y^*) = 1 + e. 

In order to the determine the dimensions of operators at V2 , we will need to know the 
functional form of 7x- Expanding jx in powers of A and yi, we find 

7x = tt 2 A 2 - 4tt 4 ^I(i/ i5 e)c iyi \ 2 + ... . (2.6) 

i 

We have relegated the detailed computations to appendix A. Here q is the OPE coefficient 
of OO^ —¥ Li in ( |2.1|) . X is a function of Vi and e which is smooth everywhere, and which 
we have normalized so that: 

lim J(i/ i} e) = 1 . (2.7) 

The ... in (|2.6| ) includes important corrections at 0(X 4 ). Ideally we would also compute 
these corrections for completeness. However, the calculation is third order in conformal 
perturbation theory and beyond the scope of this work. 



Finally, we also need the leading-order beta functions for the couplings yi. There is a 
classical contribution to the beta functions from the explicit scale dependence, j3 yi = z^j/j + 
.... The leading order correction in A was already determined using conformal perturbation 
theory in |15| ; we have redone it in appendix A with the notation and conventions of this 



paper for completeness sake. As a result of this calculation, one finds 

P Vi = ViVi ~ -^CiX 2 + . . . , (2.<< 



where Ci is again the same OPE coefficient as appears in (|2.1] ) and ( |2.6|) . We see that yi 
will flow to a fixed point y^ ~ \^-\\. 

For the theory to be under control and conformal perturbation theory to be valid, we 
obviously need y^ <C 1. This places a restriction on the range of U{ relative to the other 
parameters. In general, we need 

u { > Cl e . (2.9) 

If Ci ~ 0(1), then V{ must be parametrically larger than e. If, on the other hand, q <C 1 
(as in the perturbative examples we will study in section 3), then z/j can be much smaller. 
Of course, theories where Vi violates ( |2.9|) are perfectly valid, they just cannot be described 
using conformal perturbation theory around yi = 0. Instead, it makes more sense to treat 
the corresponding operators Li as nearly-conserved currents. In appendix B, we outline 
how to treat exactly conserved currents; we expect that nearly-conserved currents behave 
in a similar way. We will leave a complete analysis to future work. 

In summary, the RG flow of ( |2~3| ) is described by the system of beta functions 



(2.10) 



f3 x = -e\ + X ( tt 2 A 2 - 47r 4 ^ Q X(^,e) 2 / l A 2 + . . . j 



1 2 

p yi = ViVi - -CiX + ... . 
The fixed point V2 is attractive and occurs when the couplings are 

- ° l JL 

y%* — r, 9 t-\r . . . , 



t 2 e 2 ( 2 - U ) 

A 2 = — + — ^c 2 X(^,e)- + ... . 



71^ 7T^ '-r* Vi 



This is sufficient information to determine the dimensions of the non-chiral, scalar operators 
at V2, which we will turn to in the next subsection. 
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2.3. Operator dimensions 

Given the beta functions that describe an RG flow, it is straightforward to determine 
the set of scaling operators and their dimensions that describe the infinitesimal deforma- 
tions away from a fixed point. Specifically, if we deform any fixed point by J d 4 x k A 4_a E 
where £ is some operator of dimension A, the beta function for the coupling k must take 
the form j3 K = (A — 4)k + 0(k 2 ). Reversing this logic, if we are given the beta functions 
for some couplings Ki near a fixed point, we can determine the anomalous dimensions of 
operators at that fixed point by Taylor expanding and diagonalizing the beta functions 
around it, 



P* 



dKj 



( Kj - Kji! ) + 0((k - k*) z ) . (2.12) 



The eigenvectors and eigenvalues of d K f3 K . determine the scaling operators and anomalous 
dimensions, respectively. 

Applying this to the beta functions ( |2.10|) at the fixed point (|2.11| ), we find the matrix 
of anomalous dimensions at V2-& 

F = ^{ W ,A}^{ W ,A}| to =»,*,A=A. = (_ 47r 4 c ^. e)A 3 26 * J ■ (2 - 13) 

As discussed in the introduction, we are interested in the anomalous dimensions relative 
to 2'jx = 2e. It is convenient to introduce the shifted matrix 

Ar^r-2el=f h~ 2 r €) \s ~°' K ). (2.14) 

The eigenvalues 5 of this matrix satisfy the characteristic equation 

^ + 4e2 E ^! ( 2e-g ) x 11(^-26-0=0. (2.15) 



where we have used ( |2.11| ) to replace A* with \ftJT\ (again, dropping higher order correc- 



tions for now). The minimum eigenvalue 5 m i n = A m i n — 2 Ax is the quantity of interest, 
introduced in (|1.2|). In order to understand the behavior of S m i n , it will now be useful to 
distinguish different scenarios for the v%\ 



To keep the presentation from becoming too cluttered, we will simply drop the higher-order 
terms in conformal perturbation theory in the following formulas. We will discuss some of these 
corrections shortly. In the next subsection, we will analyze their effects more systematically. 



Scenario 1: i/j<1 for some i. This scenario applies to all fully perturbative SCFTs (generalized 
Banks-Zaks models), because in these theories the kinetic terms are always approxi- 
mately dimension 2. In the Vi^\ regime, two things happen. First, Z(z^,e) — > 1 as 
shown in appendix A. Then Ar is related by a similarity transform to a symmetric 
matrix 

^-=( ( ^ 4b IT) < 216 > 

The minimum eigenvalue of any symmetric matrix is always smaller than all the 
diagonal elements, so we conclude that 

Smin < (2.17) 

in this regime. The second thing that happens in this regime is that the mixing 
with X^X is enhanced, making its contribution to the anomalous dimension reliable 
compared to the 0{e 2 ) from higher order in conformal perturbation theory. The 
situation is clearest if e <C Vi <C 1 , where we find from ( |2.15| ) 



4e 2 c 2 

Srmn « L • (2.18) 

This is indeed enhanced relative to the 0(e 2 ) corrections to the beta functions which 
we have not computed. If on the other hand e ~ j/j < 1, then ( |2.18| ) no longer applies 
and we have to do degenerate perturbation theory instead. But, in this case, it is easy 
to see that the enhancement is strictly larger - instead of ( |2.18| ), 5 m i n < scales as 
(9(e). Finally, if Vi <C e, then there is already an operator Li in the OPE of X^X with 
much smaller anomalous dimension than 2e, so 5 m i n < trivially.H 

Scenario 2: V{ ~ 1 for all i. Theories in this class necessarily contain some strong dynamics. When 
all the operators Li have 0(1) anomalous dimensions, the mixing with X^X is not 

enhanced, and we find 

o *— \ c 2 X(i/,-, e) 

<W^-4e 2 V J l r } (2.19) 

j J 

from ( |2.15| ). In this case, we cannot determine the sign of 5 m i n . Even if we knew the 
sign of X(z/j, e), it turns out higher order corrections in conformal perturbation theory 



Technically, this assumes that the corrections to Ui from mixing are parametrically smaller 
than Ui itself. If instead quantum corrections to Vi are large, then conformal perturbation theory 
is clearly breaking down. We will discuss this further in the next subsection. 



are equally important to the contribution ( p.!9|) . In particular, 5 m i n is sensitive to 
(9(A 4 ) terms in (3\ which we have not computed (the ... in ( p.lOQ ) . Nevertheless, 
we see that the lowest dimension operator in the X^X OPE has dimension A m ^ n = 
2 Ax + Smin = 2 + 2e + 0(e 2 ) to leading order. This is consistent with the bounds 
of pnj] , who found that A min < 2 + 2e + 2.683e 2 . It would be very interesting to 



see whether this bound is saturated after including the (9(A) contributions. If this 
bound is the strongest possible one, then it may be possible to compute the numerical 
coefficient at 0(e 2 ) directly from conformal perturbation theory. 

To summarize: in this section, we found that the lowest dimension operator that can 
appear in the OPE of X and X^ has dimension that is at most A m j n = 2 Ax + |(9(e 2 )|. 
When Vi <C 1 for some i, as is the case for all fully perturbative SCFTs, we found that 
A m j n — 2 Ax ~ — e 2 /^i < 0. The bounds when Vi <C 1 are stronger than those found in 



TI}| . It would be interesting to determine all 0(e) corrections to see if this is true in 



general. 

2.4- Taking into account higher-order corrections 

In the previous subsections, we determined the fixed point and the matrix of anoma- 
lous dimensions by neglecting higher order corrections in conformal perturbation theory. 
Now let us discuss the conditions under which this is valid. 

One basic assumption which we have been making implicitly so far is that there 
are no large hierarchies in the OPE coefficients. Without this assumption, leading-order 
conformal perturbation theory is not necessarily valid. For instance, if the coefficient of 
y\ 2 in ( |2.6|) is anomalously small, then higher orders in conformal perturbation theory 
involving unsuppressed OPE coefficients could dominate. 

In other words, we are assuming that parametrically, ([2.6|) and (|2.8[) are really 



P Vi = "iVi + « ( -ttqA 2 + . . . 



7 x = tt 2 A 2 + K Vi (-47r 4 X(^, e)5,A 2 + ...)+ k 2 \\. . .) 

1_, \ (2-20) 

— < 

2 

Here Ci and all the coefficients in . . . are assumed to be 0(1). k can either be (at most) 
(9(1) or parametrically small. In the BZ models of Section 3, k ~ 1/N. 

Having clarified this point, let us now systematically consider higher-order corrections 
to 7x and j3 Vi , the ... in ( |2.6| ) and (|2.8|). Corrections to 7x must be proportional to A 2 ; 
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combining this with ( |2.20|) , we conclude that they are always subleading for determining the 
fixed point couplings. As discussed in the previous subsection, they are also subleading 
for determining the anomalous dimensions, provided that there are L, operators with 
dimensions close to two. All higher-order corrections to f3 y . proportional to A 2 are also 
subleading for determining the fixed point couplings and the anomalous conditions, in 
general. 

This leaves corrections to (3 y . which are O(X ), i.e. which involve the y's only. Since 
these should be compared with the classical term i^yj, they can conceivably be important 
when Vi <C 1. For instance, consider an 0(y 2 ) correction to (3 yi a 

S 0Vi = ^Cijkyjyk (2.21) 

These will affect the matrix of anomalous dimensions as 5Tij = CijkVk ~ c ' 3fcCfc<: . So if 
Vi is too small, even if it satisfies (|2.9|) , this perturbation to the matrix of anomalous 
dimensions could be larger than the classical dimension. For this paper, we will simply 
assume that these higher order terms are subleading. Nevertheless, it is plausible that 
this assumption could be promoted to a consequence of y^ <C 1. Because the two-point 
function ((Q 2 Li)(Q 2 Li)) oc V{, we must have ((Q 2 Li)(Q 2 Lj)0) oc ViVj in order for the 
(Q 2 Lij)0 OPE coefficient to be finite as Vi — >■ 0. The resulting v(Uj suppression of the 
(Q 2 Li)(Q 2 Lj) OPE ensures that STij is negligible when yij <C 1. However, a careful 
derivation valid to all orders in perturbation theory is beyond the scope of this work. 

More generally, as mentioned already in the previous subsection, in order for our anal- 
ysis in this paper to be valid, we must avoid extremely small values of z/^. The classical 
approximation of (3 Vi = Viyi + . . . should be a reliable starting point for conformal per- 
turbation theory. For this to be true, U{ must be parametrically larger than any of the 
corrections to the dimension of Li coming from conformal perturbation theory. In the BZ 
models with k ~ 1/N being the smallest parameter, this is automatically guaranteed. 

Overall, the complication presented in this general analysis is that, a priori, all the 
OPE coefficients and dimensions appear to be independent parameters. If one of these 
numbers is anomalously small or large, it could invalidate our leading order calculation. 
In practice, a given theory (or class of theories) is controlled by a much smaller set of 
parameters. Using these parameters, it is easier to establish the regime of validity of 
conformal perturbation theory. 



These corrections seemingly could lead to new fixed points with y ~ v. However, it is easy 
to check that such fixed points are repulsive, not attractive. Around the attractive point we 
considered above, 0(y ) corrections are clearly subleading as far as determining the value of y* 
is concerned. 
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3. Banks-Zaks Fixed Points 

In the previous section we presented general results on the spectrum of operators 
in the OPEs of a nearly free, chiral superfield X. Now, we will illustrate these general 
results with broad class of simple, calculable examples, namely supersymmetric Banks- 



Zaks |35|,[| (BZ) fixed points and their generalizations. We will examine the dimensions 
of operators at BZ fixed points from a more conventional perspectives and show that the 
results agree with the analysis using conformal perturbation theory. 

Before we start with a detailed analysis, let us make a few general remarks about 
perturbative SCFTs. In four dimensions, gauge fields are required in order to produce a 
non-trivial, perturbative fixed point. This requirement has several important consequences: 
(1) Charged matter cannot produce chiral operators of dimension near one because the 
matter fields are not gauge invariant. (2) To produce an operator of approximate dimension 
one, we must introduce a singlet, X, under the gauge group; therefore, the only marginal 
or relevant interactions are through operators in the superpotential. (3) The kinetic terms 
for the matter fields are approximately dimension two and will appear in the OPE of X 
and X\ (4) Therefore, all the models of this type are described by "scenario 1" of the 
previous section and we must have S m i n < on general grounds. The two examples in this 
section will serve as illustrations of these basic observations. 

3.1. Simple Banks-Zaks fixed points 

The simplest supersymmetric Banks-Zaks fixed points arise in M = 1, SU(N C ) gauge 

theories with Nf = 3N C /(1 + e) flavors (we assume Nf, iV c — > oo throughout). The 

only coupling in these theories is the gauge coupling; since we are at large N, we will use 

instead the 't Hooft coupling 

._N e g 2 

9 = ^- (3 - 1} 

Its RG flow is described (in a suitable scheme) by the NSVZ beta function ||37|| : 

^3 (3JVc-JV>(1-27q) 



& 



->• 



167T 2 1 - iVV/87r 2 (3 2) 

_6 f_ 

1+el-g 



^ = -1^1^(1+70)— 6WCd) (| +7fl) 



In particular, we will take formulas and conventions for perturbative beta functions and 



anomalous dimensions directly from [36]. 
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where we have defined f(g) so that f(g) = 1 + . . ., and 

7Q=("^ + ^ 2 )) + ^(^ + ^ 2 ))+^(^) (3-3) 

is the anomalous dimension of the fundamental flavors. We have written the anomalous 
dimension in the form ( |3.3| ) to emphasize the nature of the double expansion in 1/N C and 
g. The Banks-Zaks fixed point occurs at j3§ = 0, which according to (|3.2|) is equivalent to 

1Q = -~e , (3-4) 

or in terms of the 't Hooft coupling, 

9:0 = (c + Otfj) +^(e + O(e')) +o(±). (3.5) 

Here the subscript on g* is to emphasize that this is the fixed point in the undeformed 
BZ theory. Shortly we will deform the theory by coupling it to X , and that will change 
the value of g* . 

The BZ fixed point will play the role of V\ in Section 2. We would now like to identify 
some linear combination of the gauge invariant operators QfQg with the operator O in 
section 2. According to ( |3.4| ), these have dimension close to 2: 

A(Q f Q § ) = 2-€. (3.6) 

Note that since Q and Q are not gauge invariant operators, there are no chiral primaries 
with dimension near 1. This is a general feature of all BZ-type theories. 

As discussed in Section 2, the principal requirement on the operator O (other than 
having dimension close to 2) is that after deforming by J d 2 9\XO, the U(l)x symmetry 
is completely broken, so that no conserved currents appear in the X^X OPE.H Modulo 
SU(Nf)L x SU(Nf)n rotations, this condition is uniquely satisfied with the choice 

= a5 ff ~Q f Q f . (3.7) 



Note that in Section 2, we assumed for simplicity that O was a singlet under the global sym- 
metries of the SCFT, whereas here QfQg transforms as (Q, □) under the SU{Nf)r, x SU{Nf)R 
global symmetry. This complication is dealt with in Appendix B - the upshot being that in gener- 
alized BZ theories, provided that no conserved current appears in the X'X OPE, the generators 
that are broken by the deformation do not mix with X^ X and Li at the order we study and 
decouple from the problem. This is doubly guaranteed in BZ theories, by the fact that the Li are 
singlets with respect to the global symmetry, and by the fact that we are at large N. 
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The normalization a = , """ (1 + 0(g 2 )) is chosen so that O is CFT-canonically nor- 
malized (with Q, Q canonically normalized). The deformation breaks the flavor group 
U(l)x x SU(Nf) x SU(Nf) -> SU(Nf) d i ag , so there is no unbroken £7(1) which can ap- 
pear in the OPE of X and X*. 

To make contact with Section 2, we must also determine the OPE coefficients and 
dimensions at V\. Given flT7|), the 0+0 OPE is 

O(x) f O(0) ~ |x|- 2Ao + c|a;| Ai - 2Ao L(0) + . . . . (3.8) 

Here c is the OPE coefficient, and L is given by 

L = 6Tr(Q t Q + Q t Q). (3.9) 

To leading order in perturbation theory (i.e. with free- field contractions), the coefficients 
b and c are given by 

4tt 2 / 2 , s 

6= . c=J— — . (3.10) 

y/2NJK y N f N c 

The additional scalar operators in ( |3.8| ) can be safely ignored. Any Lj must be a 
SU(Nf)diag x U(1)r singlet and be invariant under Q -H- Q, and ( |3.9| ) is the unique such 
operator with approximate dimension 2. Note that the OPE coefficient is 1/N suppressed. 
This is a general feature of BZ-type theories and is important for the decoupling of global 
symmetry currents explicitly broken by the deformation of V\. We also note that L is 
a singlet under the full SU(Nf) x SU(Nf) which leads to a further decoupling of these 
currents. 

We see from (|3.9|) that L is nothing but the anomalous Konishi current of the SQCD 



theory, which satisfies the anomaly equation D 2 L oc TrVF a VF a ||38|| . Therefore the anoma- 
lous dimensions of L and Tr W 2 are the same. The latter can be derived by expanding 
the beta function for the gauge coupling around the fixed point (see the related discussion 
at the beginning of Section 2.3). We conclude that 

VL = dgPg\g=g,, = -6^ 2 , /(^,o)^7q(^, ) 

= ^ + o^)) + l l ^ + o^)) + o( y ±) , (3 ' n) 

where f(g) was defined in (|37^), and we have used (|373|) and ( |3.5| ) in the last equality. The 
anomalous dimension of L is 0(e 2 ). 
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Finally, let us deform the theory by 

W 



A-xo = £«"«,«, 



(3.12) 



and work out the matrix of anomalous dimensions in this simple example. Here a was 
defined below ( |3.7| ). We will first compute the anomalous dimensions using a more con- 
ventional perturbative approach. Then we will compare this with the general conformal 
perturbation theory of Section 2. 

In the deformed theory, the anomalous dimensions of X and Q are given by the 
following: 



1 



1Q 



IX = 



1 



"2 I 1 N*) S+ 2N? 



1 



* + °l*'.&.£ 



3 ~ f gX_ A^ 

2(1 + e) + \N?'N? 



where we have introduced a large iV coupling for A, A 
consist of ( |3.2j) , together with 

/ % = 2A(2 7Q + 7x ) . 

While these equations are exactly solved by 



Ni a 2 X 2 



(3.13) 



The beta functions 



(3.14) 



7Q = ~2> 7x = e 



the couplings at the fixed point are determined perturbatively using ( |3.13| ), 



9* — 9*,o + 



A, 



3iV 2 



2 / e 2 



+ o 



\m) 



(3.15) 



(3.16) 



At the new fixed point, we note that g* differs from g*^ at 0(e/iV 2 ). 

Given the couplings, we find the matrix of anomalous dimensions is given by: 



9 dX 

dg d\ 



I ■ _ I dg Q\ 



g=g*,X=X* 



-6g*(9*) d g'fQ(9*, A*) -6g*f(g*)d' X iQ(g*, K) 

2Kdg(2~(Q + 7x)(^*, A*) 2A*^(27 Q + 7x)(<7*, A, 



(l + 3%) *l + Otf /N*) -fj + O(eyN^) 
- 1 6 + b(e 2 ) + (9(e/iV 2 ) 2e + ^ + O^jNl) 



(3.17) 
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The eigenvectors of this matrix correspond to scaling operators formed from L and X^X, 
with anomalous dimensions given by the eigenvalues: 

2e+! i^r +0 (w)- -+%"+°(^)-- < 3 - 18 > 

In ( |3.17| ) and ( |3.18| ), we have explicitly included all the terms that are reliably computed 
with one-loop anomalous dimensions. At the same time, we should emphasize that the i>l 
appearing in ( |3.17| ) and (|3.18| ) is the exact anomalous dimension of L in the undeformed 
BZ theory, i.e. including all higher loop corrections. It is important to capture these 
corrections even though they cannot be computed explicitly, so that the deformation only 
shifts v L by 0(l/iV 2 ). 

We can arrive at the same results using the general formalism from section 2. Accord- 



ing to ( |2.16| ), we should have 

v L 2c L e 



2c L e 2, '+- < 3 ' 19 > 



so the eigenvalues are: 



4e 2 c| 4e 2 c| 



2e+- ^ + ..., u L -- ^ + ... (3.20) 

2e - v L 2e - v L 



If we plug in cl = \ N 2 N — \ Lv 2 ' then we find perfect agreement with the correction 
to the 2e eigenvalue in ( |3.18| ). On the other hand, our 1/iV 2 correction to the vl eigenvalue 
does not match. This is because we are keeping enough terms in conformal perturbation 
theory to determine the correction to Ajtx, but not enough terms for the correction to 
A^. For instance, a term ~ c^A 2 in (3 Vi would correct vl at the same order as what one 
gets from mixing (but it doesn't affect the correction to 2e to leading order). While it 
would be interesting to also get these terms right, it is not necessary for our purposes and 
is beyond the scope of this paper. 

Of course, given (|3.11|) , we know that this example does not come close to testing 
whether the dimension of X^X can be larger than 2 Ax; in terms of the classification of 
section 2.3, we are in the v <C e regime of "scenario 1." This occurs in most BZ-like fixed 
points - the X^X OPE always contains an operator associated with an anomalous U(l), 
like L, whose anomalous dimension v^ typically obeys v^ ^C e. As we will show in the 
next subsection, it is possible to construct a model with all ul > 2e, but at the cost of 
significant tuning of parameters. 
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N 4 



Figure 1: Quiver diagram for a viable generalized BZ model, where it is possible to achieve 
ul 3> e for all L. Circles denote SU(N) gauge groups; squares denote SU(N)l X SU(N)r global 
symmetries. Links between nodes denote vector-like bifundamental fields. Finally, every triangle 
corresponds to a gauge and global-symmetry invariant Yukawa coupling. 



3.2. Generalized Banks- Zaks Fixed Points 

Clearly, to get a meaningful example with all vl > 2e, we need to consider generalized 
BZ fixed points which have multiple small parameters and the freedom to dial e — > 
without simultaneously dialing any of the v^ — > 0. 

Shown in fig. 1 is a quiver diagram for a viable generalized BZ model. Although the 
model is quite complicated, it is the simplest example we have been able to find. The 
nodes in fig. 1 are labelled by iVj, i — 1, . . . , 6; circles indicate SU(Ni) gauge groups, and 
boxes indicate SU^Ni)^ x SU(Ni)ji flavor symmetries. Links between nodes (there are 10 
altogether) denote vector-like bifundamental fields; we will label these fields by the pair of 
numbers corresponding to the nodes that they connect, e.g. the bifundamentals between 
nodes iVj and N 2 will be called Q12 and Q2i- Our convention will be that the field is 
always a fundamental under the first node and an anti-fundamental under the second 
node. Finally, every triangle corresponds to two gauge and global- symmetry invariant 
Yukawa couplings with different orientations. For instance, the 1-2-3 triangle corresponds 
to 5W = Ai 23 Trg 12 Q 23 Q3i + A' 123 Tr Q21Q32Q13. 
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The non-anomalous, unbroken flavor symmetry of this model is Yli=4^^(-^i)L x 
SU(Ni)n, together with 5 baryonic symmetries left unbroken by the Yukawa couplings. In 
addition, at the fixed point there is an unbroken Z2 symmetry under which Qij -h- Qji. 
We will choose to preserve this Z2 symmetry along the entire flow, to avoid unnecessary 
clutter. As a result, 7^ = 7^ and Xijk = \kji- One can also include the Z2 odd couplings 
and verify this explicitly. 

As we will show, in this model e = —2712 can be dialed to zero while keeping all the 
gauge and Yukawa couplings (and hence all the v£) nonzero. Therefore, we will couple X 
to 

0~TrQ 12 Q 21 . (3.21) 

Before we launch into a more detailed analysis, let us discuss some general points: 

1. O appears to be the unique operator in this model whose anomalous dimension can 
be dialed to zero independently of the couplings. 

2. Unlike in the previous example, O is a singlet under the gauge and global symmetry 
of the quiver theory. The analysis of section 2 will apply directly, without having to 
deal with any of the complications discussed in appendix B. 

3. In general, there are a number of approximately dimension two scalar operators Lj 
that could appear in the OPE of X^ and X. Clearly, any such operator must be a 
gauge and global singlet, and they must be invariant under Qij -H- Qji. The complete 
list of compatible dimension two operators is 

LlJ = jm?r Tr {Q ^ QlJ + Q ^ Qji) • (3 ' 22) 

There are 10 such operators, and they are all explicitly broken by Konishi anomalies 
and the Yukawa couplings. There are an additional ten, Z 2 -odd operators of dimension 
near two, of which five are the conserved baryonic currents. These twenty operators 
can be combined to form the twenty U(l) currents of the free theory that commute 
with the non-abelian flavor symmetry {Qij — > e ij Qij). 

4. For the choice of ( |3.21| ), only one operator, Li 2 , appears in the X^X OPE to leading 
order. Its leading order OPE coefficient is: 



c -te' (3 ' 23) 

5. This model has both Yukawa couplings and multiple gauge groups. Both are necessary 
conditions - otherwise sending 7^ — >■ necessarily forces vj_, — > faster for some L. 
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Now let us embark on a more detailed analysis of this model. We will follow the 
notation and conventions established in the previous subsection, with some minor addi- 
tions. First, because of the complexity of this model, we will not be as careful about the 
double expansion in iV and e as we were previously; we will stick to leading-order one-loop 
anomalous dimensions at large N. Second, since here we have many Ni — y oo, we will 
introduce an auxiliary parameter N, and define 

Ni h A Ngf - NX% k 

^ = lv' €l = iV' *=8^"' Xijk= ^- (3 " 24) 

where b{ = 3Ni — £\. Nj are the coefficients of the one- loop beta functions. We will take 
N — y oo holding fixed x\ and e^. 

For a general quiver of this type, the gauge beta functions are (we neglect the denom- 
inator in the NSVZ formula, because it will not matter for our discussion): 

lij ~ r > X i9i y X j9j + rf / j X k\jk ■ 

k 

Here gi = if % is not a gauge node, and \\jh = if (ijk) is not a triangle. One can check 
that all of the anomalous dimensions 7^ (and hence all of the couplings gi, Xijk) are fully 
determined by setting the beta functions to zero. 

In ( |3.25| ), the gauge and Yukawa interactions contribute with opposite signs to 7^. 
This is ultimately what allows us to tune 712 —y while keeping all the couplings nonzero 
(and real and positive). To demonstrate this explicitly, consider taking 

7 3 e 4 9 9 

^1 = -r + TTo •■£> x 2 = 2 + -£, x 3 = - + -f, x 4 = x 5 = x 6 = l. (3.26) 

4 2 2— e 5 4 5 

with e, ( C 1. The origin of coupling-space corresponds to e = £ = 0. Solving ( |3.25|) , one 
finds 

712 = ~\e (3.27) 

independent of £. Thus we can dial 712 — > while maintaining £, and hence the couplings, 
nonzero (and real and positive). Correspondingly, at e = 0, the dimensions of the L^ 
operators are all nonzero: 

A L « (18.99f , 7.57£, 6.38£, 2.78£, 2.18f , 0.37£, 8.10£ 2 , 2.25£ 2 , 0.08£ 2 ) + . . . . (3.28) 
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with . . . denoting higher powers of £. 

To complete the construction of the model, we couple X to a Tr Q12Q21 via a Yukawa 
coupling A. This introduces the beta function for A = 8 ° 2 and modifies the anomalous 
dimension of Q12, as in the previous subsection. Along with (|3.25| ), the RG flow is now 

described by 

/ % = 2A(2 7l2 + 7x ) 

ix = ^x x x 2 X (3.29) 

1 X 

7l2^7l2 + 7^A. 

At the new fixed point, we again have 7x = —2712 = e. 

So as we dial e — > while keeping the couplings nonzero, we will reach a regime where 
7x = e ^C z/^ <C 1. This in turn provides a concrete example of the nontrivial part 
of "scenario 1" from section 2.3, where mixing with the L operators gives the dominant 
contribution to 5 m i n and forces it to be negative. We can see this explicitly by using the 
formulas ( |3.29| ) to calculate the correction to the X'X anomalous dimension in the limit 
( ^26|) . For e<£ 2 , we find: 

0.57e 2 

0~min = A min — 2Ax ~ — 7^2,^2 - (3.30) 

This is the parametric behavior of 5 m i n ~ — c 2 e 2 /^ expected from Section 2, with the 
smallest u^ ~ ^ 2 dominating from (|3.28| ), and Ci ~ 1/N according to ( p.23|) . We will leave 
a more detailed matching with conformal perturbation theory to Appendix C. 

In summary, we have constructed an example that realizes all the regimes of "scenario 
1" of our conformal perturbation theory calculation. This example serves two purposes. 
First, it demonstrates that there is no fundamental obstacle to varying e and v% inde- 
pendently. Second, it illustrates how the bound derived in Section 2 arises in an actual 
example. 

3.3. Comments on double-trace operators 

In our treatment of the generalized BZ theories, we have neglected an important 
subtlety: the presence of double-trace operators. In large A^ theories, the OPE of single 
trace operators will contain double trace operators with 0(1) OPE coefficients. 

When all the OPE coefficients are comparable in size, we are justified in focusing on 
the operators of the lowest dimension, as they will have the largest mixing with X*X. 
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In the weakly-coupled BZ theories, we determined the dimension of X^X by isolating 
the mixing with single trace operators of approximate dimension 2, such as the Konishi 
currents. However, the OPE coefficients that control the mixing with these operators are 
suppressed by 0(1/ N). Therefore, when N — )■ oo, single trace operators do not mix with 
X^X and do not alter its scaling dimension. One may worry that the dominant effect 
arises from the double-trace operator L4 = O^O which has approximate dimension 4 but 
an unsuppressed OPE coefficient with O. Since 1/N is the smallest parameter in the 
model, if L4 has any unsuppressed mixing with X'X , then it will dominate the anomalous 
dimension calculation, despite being of much higher dimension 

Fortunately, the double-trace operator is also decoupled from X^X to leading-order 
in 1/N. We see this from the calculation in appendix A: the y\ 2 term in ■yx vanishes when 
Al = 2A©. Deviations from this are proportional to A^ — 2Ae> ~ 1/N 2 for the double- 
trace operator in question. The double-trace operator then contributes at the same order 
as any other real scalar operator with v ~ 0(1), and its mixing with X^X is subleading 
relative the Konishi currents which have !/<l. 

What about at higher orders in conformal perturbation theory? These correspond to 
higher order (four points and higher) correlation functions involving (9, 0t, and L4. In 
the N — 7- 00 limit, we may use L4 = O^O to find that non- vanishing correlation functions 
are disconnected and do not require additional counterterms. New counter terms arise 
from connected correlation functions of O and O^ and are necessarily 1/N suppressed, as 
standard 1/N power-counting shows. We conclude that the decoupling of the double-trace 
operators is robust to all orders in perturbation theory. 

4. Conclusions 

A very interesting question, with numerous potential phenomenological and formal 
applications, is: what are the bounds on the lowest dimension operator that can appear 
in the OPE of scalar primaries? To date, most work on this question has focused on the 
constraints of unitarity and crossing symmetry |6|-|i~l|, or on theories with gravity duals 
|3t| . In this paper, we address how such bounds arise in perturbative SCFTs. Using 
conformal perturbation theory, we found that a bound A x tx < 2 Ax can be seen directly 
from the beta functions. 

The calculations in this paper were performed as a perturbative expansion in e = 
Ax — 1. In perturbative theories, the 0(e 2 ) contribution to A x tx was computable because 

20 



it is dominated by mixing with other operators with dimensions near two. These results 
are sufficient to cover a very broad class of perturbative SCFTs, including Banks-Zaks 
theories and their generalizations. 

Extending our results to all SCFTs would require a more comprehensive treatment of 
conformal perturbation theory in several regimes: 

1. We did not include all C(e 2 ) contributions to the beta functions, including an C(|A| 4 ) 
contribution to the anomalous dimension of X . These corrections will be relevant to 
strongly coupled theories. 

2. When the deformation breaks a global symmetry, the RG flow is more complicated 
and a complete analysis might be insightful. In the appendix, we simply showed that 
the contributions to A^tx from the broken currents are small in the cases of interest. 

3. We did not completely elucidate the limit where the dimension of a non-chiral oper- 
ator approaches two. As these operators become conserved currents in this limit, a 
description using approximately conserved currents would likely be more useful. 

4. Finally, as with all perturbative arguments, our results do not apply if contributions 
at a given order in perturbation theory are anomalously small. A more systematic 
understanding of conformal perturbation theory would be helpful to address these and 
other questions. 

The calculations in this paper could be useful for understanding RG flows in other 
applications. One advantage of our approach is that we have made all properties of the RG 
flow manifest. The results in the paper are organized explicitly as a perturbative expansion 
in e and are non-perturbative in any other small parameters present in explicit examples. 
For this reason, our analysis could be applied to strongly coupled SCFTs that are weakly 
coupled to a free field. For example, many theories contain a regime of parameters where 
the dimension of a gauge invariant chiral operator violates the unitarity bound and is 
believed to become a free field |39| , |40[1 . It is plausible that our results could describe these 



theories as the dimension of this operator approaches one. 
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Appendix A. Explicit Calculation of the Beta Functions 

A.l. Setting up the calculation 

In this appendix, we will explicitly calculate the contributions to the beta functions 
quoted (2Jj) and ( |2.8|) . We will largely follow the approach taken in [15 . However, we do 
offer an improved understanding of total derivative counterterms, and to our knowledge 
we perform the first "two-loop" calculation in 4d conformal perturbation theory. 

We would like to understand the renormalization group flow of the Langrangian ( |2.3| ) . 
For convenience we will write the action as S = S-p 1 + 5S with 



+ ( f d 4 z + d 2 6^-\A e OX(z+,9) + I d A z~ d 2 ^-\A e O^X\z~ ,9) 



(A.1) 



As in the main text, we have taken Aq = 2 — e and A^. = 2 + z/j. Here we are being careful 
to show explicitly the dependence on the superspace coordinates (z^) = x M ± iOo^Q. In 
order to achieve finite results, we will impose UV and IR cutoffs A and £~ l respectively. 
We are free to choose a holomorphic renormalization scheme in which A only depends on 
the cutoff through its classical scaling, A = Ao(A/Ao) -e , where Ao is some fixed scale. 
Also, by definition, yi refers to the part of the Li coupling which has the classical scaling, 
yi=y i0 {A/Ao) u \ 

We will determine the beta functions by requiring that the correlation functions are 
independent of the cutoff scale A. The counterterms SZx and 5yi in (|A.1| ) are fixed by 
this requirement. Let us write 5Zx and 5yi as an expansion in the couplings: 

5Z X = ai|A| 2 +auyi + a 2 |A| 4 + a 2 i\X\ 2 yi + a 2 ijyiyj + ... 

(A.2) 
Syi = b u \\\ 2 + bujVj + hi\M 4 + hij\M 2 yj + hijkVjUk + ■■■ ■ 

By expanding a general correlation function in powers of A and yi and performing the OPE 
repeatedly, we can fix the counterterms order by order. 

For our calculations, we will need the 00\ XX^ two-point functions and the LOO^ 
three point function in superspace. Superconformal invariance determine these up to 
overall normalization [II]. We fix the two-point function normalization to one, and then 
the three-point function normalization is the OPE coefficient: 



(o(zf 1 e 1 )oHz^e 2 )) = ,*,_. , (x(z+,e 1 )x^z,,d 2 )) ' 



(X + )2 (2-e)' X VI. xy V2 , .y, (x + )2 
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(A.3) 



where 



{X±Y = a& - i(9 % o»9 3 - Oj^Oi) - Wi^dij = {z~y - (z+r + 7M 5 o»Qi (A.4) 



is a supertranslation invariant interval. From (|A.3| ), we can read off some superspace OPEs 
that we will need for the calculation: 



and 



x\ z ;,e 2 )x(z+, e{) -> — \— + x^x( Xl , e 1 ,e 2 ) + ... 



U{x^ 6> 3 , 6 a )0(4, e 1 ) -»- C ; 0(z+, 9,) + ... 

Li(a ? 3 > g3,g 3 )O t (^' > g2)-> , v + C ! 2+l/< Q t ^2"^2) + .- 



(A.5) 



(A.6) 



(A 23 ; 

To be precise, the OPE limit we are taking in ( |A.5| ) is x 2 — ¥ x%, 9 2 — >• 9%, 9\ — > 9 2 - Here . . . 
denote superconformal descendants, i.e. derivatives in position and in superspace acting 
on the superfield primary. These will not matter for any of our beta function calculations 
in the next subsection. 

Finally, one comment about our regularization procedure. As mentioned above, we 
will impose UV and IR cutoffs of 1/A and £ respectively. These will be imposed as a hard 
point-splitting cutoff between any two operators. It is crucial for our calculations that this 
point-splitting cutoff be implemented in a supersymmetric way. We have found that the 
easiest way to do this is to work in superspace instead of components, and to be careful 
to always impose the cutoff on supertranslation invariant intervals Xf- defined in ( |A.4j) . 
Other ways of imposing the cutoff, for instance on position differences Xij, appear to break 
supersymmetry, leading to non-supersymmetric counterterms. 

A. 2. Calculating the counterterms 

Now we are finally ready to calculate the counterterms. Let us imagine expanding a 
general correlation function in powers of yi and A. 

1. At 0(y), we find: 

. f dtxdtd^auViXiX + {y i + bu j yj)A- Vi L i ^ . (A.7) 
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The i/i term is cutoff independent, and the other terms must be absent, since they are 
cutoff-dependent. So we conclude that 

an = buj = . (A.8) 

Note that we are not allowed to use the equations of motion of the free theory to set X^X 
to zero, since X'X appears in the correlation function as an integrated operator, and the 
equations of motion can fail due to contact terms when operators coincide. 

2. At (9(|A| 2 ), we find (now and henceforth dropping the (•••}): 



(A.9) 



/ d 4 x d 4 9 (-^ai|A| 2 X t X(x, 0, 9) + b u \X\ 2 h~ Pi Li{x, 0, 9) 

+ -^|A| 2 A 2e [ d 4 zf d 4 zZd 2 9 1 d 2 9 2 OX(z+9 1 )O^X\z2,9 2 ) 
A-n z J - 

The a\ and b\% counterterms are required in order to cancel off the cutoff dependence in 
the third term when the integrated operators fuse to give X^X and Li. To see this, we 
use the OPE ( A.5 ) on the third term: 



/ d 4 z+ d 4 Z2 ~ d 2 9 1 d 2 9 2 OX(z+, 9 1 )O^X\z 2 , 9 2 ) = 

[ d 4 z+ d 4 Z2 ~ d 2 9, d 2 9 2 ( -—J-—- xiX{x 1 , 0i, 9 2 ) + fv+ ,l_ 2e _„ cM*i, Ox A) + ■ ■ ■ 

(A.10) 
Then we are free to change integration variables z^ — > -X"^ an< ^ z i ~ * x i with trivial 
Jacobian, so that ( |A.10|) becomes 



/ d 4 X+ ^^ xj d*x dH X^X(x, 9, 9)+ J d 4 X+ ± C l 2e _ ui xj d A x d*9 L(x, 9, 9) 

(A.ll) 
Performing the X^ integrals with a short-distance cutoff 1/A, plugging this back into 
( |A.9| ), and requiring that the counterterms in ( |A.9| ) cancel off the result, we finally obtain 



e 2(2e + z^J 

The appearance of the e _1 in the counter-terms here is no different than the (4 — <i) _1 
that appears in counter-terms in dimensional regularization. Although the counter-term 
will diverge in the e — > limit, the e _1 does not appear in the beta functions and RG flow 
is well behaved. 
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3. At 0(y\), we find (taking into account ( |A.8| )): 

^- yi XA e - Vi [d 4 x 3 d 4 9 3 d 4 z+d 2 9 1 L l (x 3 ,9 3 J 3 )OX(z+ 1 9 1 ) + c.c. (A.13) 

2n J 

These terms must be cutoff independent, because in our holomorphic scheme, there are no 

counterterms which can cancel them. Indeed, when Li approaches OX or O'X^ , we can 

use the OPEs in ( |A.6| ). Substituting these into ( |A.13|) yields 

-^AA 6 " 1 " fd 4 x 3 d 4 9 3 d 4 z+d 2 9 1 Ci OX(z+,9 1 ) + c.c. (A.14) 

Ztv J (X 3 T 1 ) Z+Ui 

We are free to change integration variables from x 3 — > X^ and #3 — > 6*31; then (|A.14| ) 
is annihilated by the d 4 9 3 \ integral. This change of variables also ensures that all the 
superconformal descendants in the LO OPE do not contribute to renormalization of OX 
in the superpotential. 

4. At C%|A| 2 ), we find: 
1 



;y i \x\-A-'- ,/ ' l d 4 z+ d 2 e 1 d 4 zz d 2 e 2 d 4 x 3 d 4 e 3 L t (x 3l e 3 , o s )ox{4, 9 1 )o^x^(z 2 ,9 2 ) 
+ a~ 2 " J d 4 Xl d% d 4 x 2 d 4 e 2 L l ( Xll e ll di) ( ai x*x(x 2 , e 2 , e 2 ) + ^ 2 b Xl A-^u{x 2l e 2 , e 2 )) 

+ A^" 2e I d 4 xd 4 6 (a 2l X ] X(x,9,9) + ^ 2 b 2l3 A~ v ^ Lj(x,9,9)) 

(A.15) 
The first term gives rise to cutoff dependence that must be cancelled by the counterterms. 
We already fixed a\ at a lower order in perturbation theory. Clearly, a\ does not play an 
important role - L and X^X have a trivial OPE in the unperturbed theory and 01 will 
cancel the contribution to this operator from the first term. The bu and b 2 ij terms are 
relevant for the renomalization of y^k and other non-chiral operators of Vi, but these 
contributions are negligible for our present purposes. What we are interested in is the 
X^X counterterm represented by a 2 i that is required to cancel off the cutoff dependence 
when Li OO^ — V 1 in the first term of ( |A.15| ). In more detail, using the OPEs ( |A.5| ), we 
have from the first term: 

d 4 zf d 2 e 1 d 4 z^ d 2 e 2 d 4 x 3 d 4 e 3 L t {x 3 , e 3 , 9 3 )ox(zt, 9i)o^x\z 2 , e 2 ) 

-> J d 4 z+ d 2 e 1 d 2 e 2 x^x( Xl , e u e 2 ) (A.i6) 

d 4 Z2 ~ d 4 x 3 d 4 6 3 (L { (x 3 , 9 3 , 3 )O(z+, 6 1 )0\z-J 2 )) 
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Comparing with ( |A.15 ), we see that the required counterterm isQ 



a 2l = -A 2e "^ f d 4 z 2 d 4 x 3 d 4 9 3 (Li(x 3 , 9 3 , 9 3 )0(zf, 9 1 )0\z 2 ,0 2 )) (A.17) 



where the integrals should be performed at fixed z± , $i, and 6% . Using ( |A.3|) to evalu- 



ate the three-point function, and changing integration variables from (z 2 , x 3 , 9 3 ,9 3 ) — > 
(X 23 , X 31 , 03i, 9 3 i ) (with trivial Jacobian), this becomes 

a 2l = -A 2e "^ J d 4 X+ d 4 X+ d 4 9 31 
= c i (2e + v i )(2-2e-v i )A 
-l(vi,e) 



(X+ + X+ + 2 l 9 31 a9 31 )i-^(X+)^(X+)^ 

2e-Vi f d ^23 d X 31 

'., {x^ + xiy-^ix^y+Hxt) 2 ^ 

87T 4 C, 



Vi - 2e 
where we have defined the dimensionless integral function 

(u i -2e)(2e + u i )(2-2e-u i ) K2e _ Vt f d*X+d*X+ 



(A.18) 



8^ ' J {X++X+)*-*->>*{X+)*+>'i{X+)*+'>* 

(A.19) 

In general, this is a complicated function of z/j and e, because we must take care to regularize 
the integral so that \X 23 \ > 1/A, l-X^I > 1/A, and \X 23 + X 31 \ > 1/A. In two special 
cases, things simplify: 

a. If vi <C 1, then the singularities at X 23 = 0, X 31 = 0, and X 23 +X 31 = are separately 
integrable, so the integral is dominated by the common singularity at X 23 , X 31 — > 0. 
We should be able to calculate this dominant part by doing the X 23 integral without 
regard to the X 23 — and X 23 = —X 31 singularities, and then doing the X 31 integral 
with a 1/A cutoff. Doing so, we obtain X = 1 plus higher orders in Vi and e. 

b. If Aj^ = 2Ae> (i.e. v\ = 2 — 2e), as is the case for the double-trace operators discussed 
in section 3.3, then ( A.19 ) simply vanishes. Therefore, in this case, the wavefunction 



renormalization of X is independent of yi to this order. 



Actually, we are dropping here a potential contribution from J d 4 xid 4 6i\X\ 2 LiX' X x 
J d 4 X2d 4 02ViLi. This vanishes when in — >• or when Aj^ = 2Ae>, the two special cases we 
are interested in. 
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A. 3. Calculating the beta functions 

With the counterterms ( |A.2| ) in hand, it is straightforward to compute the beta func- 
tions for the couplings, keeping in mind we are interested in the physical couplings, 

\p^ = AZ" 1/2 , yf ys = y l + 8 Vi . (A.20) 

Using the counter terms and classical running of the bare couplings, we find: 

(3 XPhys =- e A^ s + 7x A p ^ s 

h s 1 ( A - 21 ) 

(3 y p h ys = v^i + 5y'i = u t y? ys - -C;|A| 2 + 0(y 2 , A 4 , y\ 2 ) 



where 7 X = - ^J^// ■ Substituting (|A^ ), (gg) , ( pTT2|) , and ( |Q5D , we obtain for the 



2 dlogA 

anomalous dimension of X: 



lx = (n 2 - 4tt 4 X(^, 2e)c t y t )\X\ 2 + 0(X\ X z y 2 ) (A.22) 

At this order in perturbation theory, we are free to substitute \ phys , yf ys for A, yi in 
( |A.21| ) and ( |A.22|) . This completes our derivation of the beta functions used in the body 



of the paper. 



Appendix B. Generalizing to the case with global symmetry 

In this appendix, we will generalize the discussion in Section 2 to situations where O 
transforms in some nontrivial representation r of the global symmetry group G of V\ . 

By definition, when O is charged, the appropriate conserved current multiplets J a 
will appear in the C-C f OPE: 

O m (x)OU0) = ^|x|- 2Ao + Y, \x\ 2 - 2A °T^ n J a + Y, \x\ Al ~ 2A ° (cOmnpMp + • • • • 

a i,p 

(B.l) 
Here m, n and p are global symmetry indices. T^ are the generators of G in the represen- 
tation r; their normalization is fixed by the requirement that the J a, s are CFT-canonically 
normalized. We are free to choose g mfl = 5 mfl at Pi. H In writing ( [B.l| ), we have allowed 



The meaning of this "metric" becomes more clear when V\ is a point on a conformal manifold. 
In that case, g m n defines a metric on that space and can have non-trivial dependence on the exactly 



marginal couplings [42]. 
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for the possibility that Li transforms in some representation h of G. Allowing for charged 
O and Li complicates the situation considerably. 

As before, we will deform the theory by the superpotential W D -^XOX, where now 
A transforms in the r representation of G (the global symmetry indices on A and O are 
implicit). This deformation breaks the global symmetry, allowing the currents to acquire 
anomalous dimensions and mix with J x = X^X and Li. 

In order to arrive at a non-trivial bound, we should forbid the appearance of conserved 
currents in the OPE of X and X< at TV That is, we will demand that our theory does 
not contain an unbroken U(l) current formed from a linear combination of J x and J a . 

These satisfy 

D 2 J a = XT a OX + ... , 

(B.2) 
D 2 J X = XOX + ... . 

Therefore XT a and A should be linearly independent vectors. Stated more formally, we 
require that the only real solution to D 2 (J2 a C a J a + CxJ X ) = is C a = Cx = 0. Using 
( |B.2| ), this requirement can rewritten (to leading order) as 

J2CaXT a + C x X = => C a = C x = 0. (B.3) 

a 

After deforming the theory by XOX, the theory will undergo an RG flow to the new 
fixed point TV As before, we will work in a holomorphic basis where the only renormal- 
ization is of operators in the Kahler potential. We should therefore consider the following 
action 



c = c Vl + fd 4 e (^z x J x + Y, z « ja + 5> 2_Ai ^ ■ L * + 



+ 



d 2 6^-A £ XOX + h.c. 

Z7T 



;b.4) 



At one loop, 5Z a = -2 7 r 2 AT a A t log(A) 0. 

The beta function for yi is calculated in the same way as before: 

Pv< = "iVi - \ Xc ^ + °( X ^ v x ^ y 2 ) ■ ( B - 5 ) 

The beta function for A is not useful in the holomorphic basis, as the fixed point appears at 
A — > oo. To study the fixed point, we will move to a non- holomorphic basis by absorbing 
Z a and Zx into the definition of A. This can be accomplished using ( p.2| ) to rewrite 

A -> A - X -Z X X - \z a XT a . (B.6) 
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The RG flow for A is governed by 

/3 A = A(-e + 7x)+7aAT a , (B.7) 

where 7 a = — | gwnf • We compute the contributions to Zx as before to find (at small 

7x = ,r 2 AAt - 47r 4 A( yi • c,)At + 0(y 2 A 2 , A 4 ) . (B.8) 

We will not compute all the contributions to 7 a , but we will simply note that it takes the 
form: 

7a = yr 2 AT a A t + 0{\\ yA 2 , y 2 ) . (B.9) 

From ( B.3|) , we see that /3\ = requires^ that 

ix = e, 7a = . (B.10) 

To leading order, this implies AT a A^ = for every a. 

We now wish to compute the scaling operators around the fixed point as before. In 
general, A is an r-dimensional vector of couplings that admits many independent deforma- 
tions. We are only interested in those "directions" that correspond irrelevant operators, 
namely the broken currents J x , J a . Other deformations either correspond to exactly 
marginal operators or total derivatives. In order to focus on just the broken currents, let 
us specialize to the deformations around the fixed point which correspond to them via 

\ = K + r] K + r] a KT a (B.ii) 

with r/o, r\ a real. In other words, we are trading the vector of couplings A for a subset 
parametrized by real (rjo,r] a ) which correspond to the broken currents. Note that the 
unbroken currents automatically drop out of this since they satisfy A*T a = 0. 

We would like to know the beta functions for (r)Q,rj a ) to linear order. We can read 
them off from ( |B.7| ), using the linear independence of A* and \*T a and the fact that 
-e + 7x, la ~ 0(rj): 

A70 = -e + ix + 0( V 2 ) = -e + 7r 2 AA t - 4tt 4 A( Ci ■ yi )X^ + 0{\\ y 2 X 2 ) + 0{r, 2 ) 

(B.12) 
Pva = 7a + 0( V 2 ) = 7 r 2 AT a A t + 0(\\ 2/A 2 , y 2 ) + 0( V 2 ) 



This conclusion is true, even when higher order terms in ( |B.2| ) are included. The solutions 
to /3\ = are equivalent to the solutions of D 2 (^2 laJ a + (7.x — e )J X ) = 0> so our conclusion 
follows from forbidding conserved currents. 
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Combining this and ( |B.5|) , we finally obtain the matrix of beta function derivatives: 



/ d PvQ d Pv 9Pr,Q \ 

dr] dr] a dy t 

dPvb d Pv b d Pv b 



9P Vl 



dVa 



\ dr)o dr) a 

2e 




v=0,yi=yi, 



dyi 
dyi / 

-47r 4 Ac i A t 

2n 2 \{T\ T a }At £>(A 2 ,y) 

Vi 



;b.i3) 



+ 0(\\y\ 2 ,y< 



■\ Ci \i £>(A 2 ) 

We can diagonalize the matrix 27r 2 A{T b , T a }\^ — >■ i> a 5 ab , and shift by 2e, to write this as 



Ar 



-47r 4 AciAt 

u a -2e 0(X 2 ,y) ] + 0{\\ yA 2 , y 2 ^ 

-Ac,At C(A 2 ) i/i - 2e 



^.14) 



In a general model, the mixing of Jx with J a is potentially important. To get a feeling for 
this, let us simplify things somewhat by ignoring the distinction between A and y and just 
taking A 2 ~ y ~ e. Also, let us drop the higher orders in ( p.!4|) , as these are irrelevant for 



our discussion. Then the characteristic equation for the eigenvalues of (p.!4| ) is 

(D a -2e-S) Uvi - 2e - 5)5 + 4yr 4 (Ac l A t ) 2 ) + C(e 2 ) x^O (B.15) 

The last term in ( |B.15| ) can be traced back to the 0(X 2 y) and 0(y 2 ) corrections to 7 a 
which we have not computed. If v a is not close to 2e, this term is higher order in e and can 
be safely ignored, meaning the leading order for 5 is as before. However, if v a = 2e + 0(e 2 ), 
then the last term in ( p.!5| ) cannot be ignored for determining the leading order 5. 

To summarize so far: we have seen that it is valid to neglect the mixing with the J a 
operators, as long as their anomalous dimensions are not extremely degenerate with 2e. If 
v a — 2e ~ 0(e 2 ), then the mixing with J a cannot be ignored, and more study (beyond the 
scope of this work) is required. 

We will now conclude our discussion by describing two important classes of models 
where the mixing can always be neglected, regardless of the J a anomalous dimensions. 

Case 1: Li are singlets under G. The correlation functions of Li with J a are restricted by 
the representation of Li under G. When all the Li are singlets of G, Ward identities 
require that (LiLjJ a ) = (or equivalently that J a is not in the OPE of Li and Lj). 
As a result, the higher order correction to 7 a take the form 



7a = 7r 2 AT a A t + C(A 4 , yAT a A f ) . 
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(B.16) 



Furthermore, yi is a flavor singlet and therefore 



^ = wi - \ c ^ + °( x ^ 2M 2 > y 2 ) ■ ( B - 17 ) 

with the Cj's now just numbers, not matrices. As a result, 

2e -4yr 4 c l AA t \ 

v a ° +0(e 2 ) . (B.18) 

-CiXX^ Ui J 

Here we have used XT a X^ ~ e 2 . Up to order e 2 corrections, u a are eigenvalues and 
do not mix with Jx and Li. So after diagonalizing the u a , v^ block as above, one 
finds that the contribution from mixing with v a is subleading in A, y, even if they are 
completely degenerate with 2e. So the results of Section 2 are unchanged. 

Case 2: Large N limit. In the Banks-Zaks examples we studied in Section 3, the OPE co- 
efficients are controlled by a separate small parameter, 1/iV. More generally, if the 
theory has a large iV parameter which controls the sizes of the O'O OPE coefficients, 

(O^OL) ~ (O^OJ) ~ 1/N (B.19) 

and such that other OPE coefficients are also suppressed by powers of N, specifically 

{LiLjJ a ) ~ l/N a , a > (B.20) 

then the conclusions of section 2 are robust. (This should include all generalized 
BZ theories and also theories with weakly-coupled gravity duals.) Using ( B.19|) and 
( BT20| ), the large iV counting of ( ETT3|) is 

/ 2e 0(1/JV) \ 

T= v a 0{l/N 1+a ) +... . (B.21) 

\0(1/N) 0(1/N) v, ) 

The 0(l/N 1+a ) came from two contributions to 7 a - the contribution at 0(yX 2 ) 
involves OO^L^ — > J a , which involves two OPE coefficients; and the contribution at 
0(y 2 ) arises from LL —¥ J, plus the fact that y ~ 0(1/ N) from (|2.11|) . So we conclude 



that the correction from mixing with v a is suppressed by l/iV 2+a and vanishes relative 
to the 0(l/iV 2 ) contributions computed in Section 2. As a result, the conclusions are 
unchanged. 
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Appendix C. More on the Generalized Banks-Zaks Example 

In this appendix, we would like to make contact between the analysis of section 2 and 
our generalized Banks-Zaks example, by rederiving ( |3.30| ) using the general formula ( |2.18|) . 

In order to use the results of section 2, we must determine the operators in the OPE 
of O oc Q12Q21 and O' along with their OPE coefficients and dimensions. This example is 
controlled by the large iV expansion and we will only be concerned with results at leading 
order in 1/N. As discussed in section 3, the only dimension approximate dimension two 
operator in the OPE is O(x)O^(0) — > L i2 where 

L 12 = ^= Tr (Q\ 2 Q 12 + QlQ 21 ) , (C.l) 

with OPE coefficient 



c - = ^- (a2) 



The problem with directly applying ( |2.18| ) is that L\ 2 is not a scaling operator in the 



original undeformed model. So instead of (|2.18| ), we should have 

5 mm = -4e 2 c? 2 (r- 1 ) 12;12 + ... . (C.3) 

where T is the matrix of anomalous dimensions in the Lij basis. The last remaining 
step is to find the transformation into this basis starting from the matrix of anomalous 
dimensions obtained from the beta functions ( |3.25|) . This change of basis is implemented by 



the Konishi anomaly and U(l) equations of motion. Adding VijLij to the Kahler potential 
is equivalent to rotating Qij and Qji by the axial phase 

Qij -> e^MQij, Qji -> e'wtQji (C.4) 

where the charge q^ oc 1/^/XiXj can be read off of the (Q^QijLij) three point function. 
This in turn induces a change in the superpotential and gauge couplings, via the current 
non-conservation equation (for a nice discussion of this, see 



The third and fourth terms correspond to the Konishi anomaly terms, and they are present 
only if the ith or jth node is gauged, respectively. Finally, from (|C.5| ), we read off the 
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shift in the couplings (our conventions for the gauge kinetic term are the standard ones, 
in which CD -±zJd 2 9TY W%) 

8£ = y^l d A 9 yjjLjj ->• SX ijk = 2(y ij q ij +y ik q ik +y jk q jk )X ijkl 5gi = -2^2y ij q ij x j g- 

i>3 J 

(C.6) 
So this is a linear map S which takes the couplings yij to the couplings (#, A). If we then 
compute T using the beta functions ( |3.25| ), we obtain it in the Lij basis by acting with 
T — y S~ 1 TS. Substituting this into ( |C.3| ), we find perfect agreement with (|3.30| ). 
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